Generation of Gravitational Radiation in Dusty Plasmas and Supernovae 
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We present a novel nonlinear mechanism for exciting a gravitational radiation pulse (or a gravi- 
tational wave) by dust magnetohydrodynamic (DMHD) waves in dusty astrophysical plasmas. We 
derive the relevant equations governing the dynamics of nonlinearly coupled DMHD waves and a 
gravitational wave (GW). The system of equations is used to investigate the generation of a GW 
by compressional Alfven waves in a type II supernova. The growth rate of our nonlinear process is 
estimated, and the results are discussed in the context of the gravitational radiation accompanying 
supernova explosions. 



PACS numbers: 04.30.Db, 97.60.Bw, 98. 62. En 

It is well known that there exist numerous mechanisms 
for the conversion betweengravitational waves (GWs) 
and electromagnetic waves [H HSBfiLaLllaLA LLil UJl 
HHQHEElElllFor example, the propaga- 
tion of GWs across an external magnetic field gives rise to 
a linear coupling to the electromagnetic field 1] , which 
may lead to the gravitational wave excitation of ordinary 
electromagnetic waves in vacuum, or of magnetohydrody- 
namic (MHD) waves in a plasma [2j,|3|,|j|. Furthermore, 
various nonlinear coupling mechanisms give rise to three- 
wave couplings between GWs and electromagnetic waves 
in matters. We also note that four-wave processes may 
cause graviton-photon conversion even in the absence of 
external matters or fields ||. Moreover, GWs can cou- 
ple to other types of waves, e.g. sound waves, also in 
neutral media [6j. There are numerous motives for con- 
sidering wave couplings involving GWs. In some cases, 
the emphasis is on the basic theory 0, 0, l^-lj^lict . In 
other works, the focus is on GW detectors [nTHaE!, 
on cosmology Hi 0, 0| , or on astrophysical applica- 
tions such as binary mergers [l7| , gamma ray bursts [lij 
or pulsars [l^. Many of the previous works have con- 
centrated on the conversion from GWs to electromag- 
netic waves, which can be analysed within a test matter 
approach which neglects the back reaction on the grav- 
itational field. We note that such an approach can be 
justified if the background energy density is low. 

In this Letter, we consider the three-wave coupling be- 
tween two DMHD waves and a GW, including the ef- 
fects of dust particles |2fjj in a dense medium such as 
the supernova where electrons, protons, and charged dust 
macroparticles are abundant. For this purpose, we derive 
the dust Hall MHD equations [2fj|. i.e. equations describ- 
ing the dust MHD waves, including the effect of a GW. 
We emphasize that for a low-beta plasma, the system 
of equations has a structure which can describe both a 
dust-dominated plasma, as well as an ordinary Hall-MHD 
plasma (if we replace the dust mass density by the ion 
mass density). Using the normal mode approach [2lj . 



the three-wave coupling equations are derived, includ- 
ing the back reaction on a GW from the Einstein equa- 
tions. The system is shown to fulfil the Manley-Rowc 
relations [2lJ (which means that the interaction process 
can be viewed quantum mechanically) and to be energy 
conserving. The three-wave equations are then used to 
analyse the generation of a GW by the compressional 
Alfven waves in the iron core of the type II supernova 
|22l Eij ] . It turns out that the characteristic timescale for 
the Alfven-GW conversion can be less than a millisecond, 
which implies that the mechanism is potentially relevant 
for the high-frequency part (> IMHz) of the supernova 
GW spectrum. 

The plasma dynamics, due to the response to a gravi- 
tational wave 

ds 2 = -dt 2 + (l + h + )dx 2 + (l - h+) dy 2 +2h x dx dy+dz 2 

(1) 

can be formulated according to 



d t n s + V • (n s v s ) = 0, 



(2) 



and 



m s n s (d t + v • V)v = -Vp s + q s n s (E + v x B) + m s n s g s , 

(3) 

where we have introduced the tetrad eo = dt, ei = 
(l-h+/2)d x h x /2d v ,e 2 = {l + h+/2)d v - 
h x /2d x ,e3 = d z , and V = (ei,e2,e3). Moreover, 

gs = - v sz ){v sx d t h + +v sy d t h x )ei 

-5(1 - v sz )(v sx d t h x - v S yd t h + )e 2 
~\ M x - v 2 y )d t h + + 2v sx v sy d t h x ] e 3 (4) 

represents the gravitational acceleration of the particle 
species s due to the GWs. We have assumed that d z w 
— d t holds for the GWs. 

The electromagnetic field is determined through the 
gravity modified Maxwell's equations. Using the same 
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notation as above, they take the form 

<9 t E = VxB-^g s ji s v s -j E , (5) 

S 

d t B = -VxE-j B! (6) 

with the constraints V • E = ^ s q s n s and V ■ B = 0. 
Here, the effects of the GWs £p) are represented by the 
effective currents 

j E = -\[{E X -B y )d t h + + {E V + B x )d t hy]e x 

-\ \-{E y + B x )3 t h + - [E x - B y )d t h x ] e 2 , (7) 

j B = -±[{E y + B x )d t h + -(E x -B y )d t h x }e 1 

-\ [(E x - B y )d t h+ + {E y + B x )d t h x ] e 2 . (8) 



and 



With the general setting established above, we will 
from now on focus on the case of a three-component 
dusty plasma, for which we have the equation of state 
p s = kBT s n s . Thus, the plasma is composed of electrons 
(e), ions (i), and dust particles (d). The mass rtid of 
the dust particles is assumed to be much larger than the 
electron and ion masses, viz. m e and rrii, respectively. 
We will assume that the plasma is approximately quasi- 
neutral, i.e. QiUi = en e — qdUd- Moreover, the waves un- 
der consideration are supposed to propagate with phase 
velocities much smaller than the speed of light c. Thus, 
we may neglect the displacement current in Ampere's law 
©, i.e. 



V x B = > a s n s v a +j E - 



(9) 



Due to the constraint m e , rrii <C md the momentum con- 
servation equation for the inertialess electrons and 
ions becomes 



= -k B T e Vn e - en e (E + v e x B) + m e n e g e , (10) 



and 



= -fcsTjVni + gin;(E + v, x B) + m^g,, (11) 

respectively. Adding Eqs. (|10f> and (|llf> . using the 
quasineutrality condition, assuming that the number 
densities of the electrons and ions are not much larger 
than the number density of the dust, and using the heavy 
dust approximation, the dust momentum equation takes 
the form 

Pd {d t + v d • V) v d = -k B (^Td - ^ T ij Vn rf 
+ (VxB)xB-jfixB + PdSd (12) 

where pd — rridUd- In Eq. l(T2^ we have also used the 
approximation [T e + (e/qi)Ti]n e < [T d - {qd/qi)Ti]n d - 



Again using Eqs. (|10|) and (|llfl to eliminate the electric 
field, Faraday's law © becomes 

d t B = V x (v d x B) - ^ V x [(dt + v d ■ V) v d - g d ] -j B , 

<ld 

(13) 

where we have used the dust momentum equation (|12fl . 

Thus, Eqs. (|T^|l and ffT^i together with the dust conti- 
nuity equation 



d t pd + V • (pdv d ) = 0, 



(14) 



constitute the dust MHD equations in the presence of a 
GW. For a low-beta plasma, the pressure term in (|12(l is 
negligible, which means that the structure of Equations 
(|12[1 - H14[) is the same as in an ordinary Hall-MHD plasma 
without dust. Henceforth, we will consider a low-beta 
plasma, drop the index d on all quantities and thus let 
q/m be either the ion charge to mass ratio or the consid- 
erably smaller charge to mass ratio of the dust particles. 
As a result, our mathematical analysis below will then 
apply either to a dust Hall-MHD plasma, or to an ordi- 
nary Hall-MHD plasma without dust. 

To simplify the problem, we consider the case when 
the dust-acoustic speed c s = fc^T/m is much smaller 

1/2 

than the dust Alfven velocity Ca — {Bq/hqp) such 
that the pressure term in ltT2*|l can be neglected. As a 
prerequisite for the nonlinear calculations, we first study 
the linear modes of the system <|12[1 - H14[) omitting the 
gravitational contributions. Letting B = B<yz, + Bi, p — 
Po + pi, where the index 1 denotes the perturbation of 
the equilibrium part, and linearizing Eqs. (|12|I - H14|I and 
Fourier analysing, we readily obtain the dispersion rela- 
tion 



k 2 z Cl) (c 2 - k 2 C 2 A ) - 



0, 



XT. 



(15) 



where lo c — qBo /m is the gyrofrequency. For frequencies 
much smaller than the gyrofrequency, we note that the 
modes separate into the shear Alfven wave, lu 2 — k 2 C\ « 
0, and the compressional Alfven wave, co 2 — k 2 C\ w 0. 
Below we will consider the more general case described 
by G3 ) however. For later applications it is convenient 
to use the linear equations to express all quantities in 
terms of a single variable. Thus, we let the wavevector 
of the dust MHD waves lie in the x — z-plane, and write 



B. t 

By 

B z 



co c ( W a - klC 2 A ) 



-iBc 



v x , v z = 0, pi = po- 



Wc {^-kic 2 A y x 



Q my 2 A x 



-(16) 
(17) 
(18) 
(19) 
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Next, we consider a system of three weakly inter- 
acting waves. Two dust MHD waves with frequencies 
and wave- numbers (wijki) and (w 2 ,k 2 ) respectively, and 
an arbitrarily polarized gravitational wave propagating 
along the z-direction with the frequency and wavenum- 
ber (ajg,k g z). Noting that the gravitational dispersion 
relation reads cu g — k g c and that Ca <C c, the frequency 
and wavenumber matchings (energy and momentum con- 
servation) can be approximated 



0«ki 



(20) 



We will thus use ki rs — k 2 below, and we define 
ki x = —k 2x = k x as well as k\ z — —k 2z = k z (let- 
ting k y = for convenience). All quantities are now 
assumed to be superpositions of two dust MHD waves 
waves whose amplitudes are weakly varying functions of 



time, i.e. p = p + Z)j=i exp[i(k 



where c.c stands for the complex conjugate, In princi- 
ple, the gravitational wave should also contribute, but we 
note that within a fluid model the gravitational wave con- 
tribution to all plasma perturbations (velocity, magnetic 
field and density) are second order in the gravitational 
wave amplitudes, provided that the GW propagates par- 
allel to B , as we have assumed. Thus, the only linear 
perturbations due to the gravitational wave are those of 
the metric as described by Eq. Q. 

Next, in order to simplify the algebra, we introduce 
the normal mode Oj defined by 



h 2 h 



B 



-D 



2p f<2 



c 2 



k]C A ) 



Bn 



-B z 



-k 2 C 



oj c k X jBo 



-B„ 



(21) 



Returning to Eqs. (|12[l - l|14|) and including the nonlin- 
ear terms 26], we can, by keeping the part varying as 
exp[i(ki • r — Wit)], derive 



dai 
~dt 



v 2x h + +v 2y h > 



uiklC\(v* 2x h y 



V *2yh+) 



2uj c (uif 



2k 2 k x B 



+ i 



k z [{u 



(B* 2x h, 



k 2 C\) 



B* 2y h y 



where 



po 
2 



k 12 C A 



J l,2 



1,4 f<2 



^21,2^1) 



V x =i 



(^1,2 

u:\k\C\k 2 — uj 2 u 2 k 2 z 



LO c (uj 2 -k 2 zl C A )k 




and 



uj x u 2 k 2 z C\ {k 2 z C\ + Wicj 2 ) u\u 2 k\ 



k 2 C 2 A )(co 2 



k z C A ) 



(27) 

Next using the Einstein equations, linearized in h + ,h x , 
keeping only the resonantly varying part of T^ v we obtain 
for the x and +- polarization, respectively 



and 



dh x 
'~dT 



IUJ 



dh y 
~dt 



P0VxlVy 2 



B X 2By\ 

Mo 



(28) 



PO {p x \Vx2 ~ V y \V y2 ) + 



B X \Bx2 — By\By 2 



/to 



which is reduced to 

dh x 
~~dT 



and 



dh + 
~dT 



W a 



V X V X \V X 2, 



V + V x lVx2, 



(29) 
(30) 

(31) 



where W g = lu 2 /2k. The total wave energy is Wtot = 



W- 



l\Vxl\ +W 2 \V X 2\" + 



W g {% 



\h + \ ) and it is eas- 



(B 2x h x - B 2y h A 



k 2 0Jiu> c k x B 

and we obtain a similar result for da 2 /dt by letting 1 <-»• 2. 
After some algebra, using Eqs. (|16|I - H19|I and J2J, we 
find that Eq. (|2*2*)l reduces to 

dv x i 



(22) 



_ ppvi r 

dt ~ W 1 '' 
and similarly for mode 2 
dv x2 p Q uj 2 



j:2 



(V+h+ + V x h x ), 



dt 



W 2 



v* xl (V + h + + V x h x ). 



ily verified from fi fy -fi ty together with (JSJ-EU that 
Wtot is conserved. Furthermore, the appearance of the 
same coupling coefficients V x , V+ in Eqs. I|23|l - I|24|) 
as well as in (|30|) - (|31ll assures that the Manley-Rowe 
relations are fulfilled, which implies that each mode 
changes energy in direct proportion to its frequency, i.e. 
(dWi/dt)/(dW 2 /dt) = uoi/u 2 etc. The system of ffity- 
({2"4"|l together with l(37)|) - ({3"l"1) describing the energy conver- 
sion between DMHD waves and GWs is one of the main 
results of the present letter. A more elaborate calcula- 
tion scheme, including effects such as inhomogeneity and 
background curvature, is a project for future research. 

We now apply our results to the gravitational radia- 
tion arising from the iron core of a type II supernova 
(23) where the densities can be of the order 10 17 kg/m 3 j2^ . 
The large neutrino outflo w ( which can reach powers of 
10 33 W/cm 2 see e.g. Ref. |22j |) can generate MHD waves 
described by (|15fl . From the flux conservation, we expect 
the iron core to be strongly magnetized (comparable to 
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pulsars), and for magnetic field strengths Bq ~ 10 s T, 
the gyrofrequency will be much larger than all other fre- 
quencies of the problem. The dispersion relation then 
separate into the shear Alfven waves ui 2 — k 2 C\ ~ and 
the compressional Alfven waves ui 2 — k 2 C\ 0. Assum- 
ing that the pump MHD wave is a compressional mode 
with a frequency 5 MHz, the matching conditions l|20() 
can be fulfilled for a GW with a typical frequency 3MHz 
and a shear Alfv en wave with the frequency 2 MHz • 
For the assumed geometry, the MHD waves couple only 
to the h x -polarization (to a good approximation), and 
by combining (|2"4"|) and i(3"0"|) . we obtain 

d 2 h* u 2 \V*\ 2 l(!,nG 2 

-&- = - h *17 a w 2 c 2 PoM ■ (32) 

Thus, noting that the factor w 2 \ V X \ /to „]¥■?. is negative 
|27| and has a magnitude of order unity [28| for the given 
parameters, the growth rate is 

7 ~yi6^^, (33) 
c 

which for a weakly relativistic pump quiver speed, 
bzil/c ~ 1/10, implies 7 ~ 10 kHz. Thus, we deduce 
that excitation of a GW by MHD waves is a reasonably 
fast process in a dense medium such as the supernova iron 
core. Furthermore, the present process can give rise to 
GWs of higher frequencies than many of the previously 
considered excitation mechanisms, see e.g. |24j . Thus, 
our model contributes to the understanding of gravita- 
tional radiation emissions accompanying supernova ex- 
plosions 



* Also at: Centre for Fundamental Physics, Ruther- 
ford Appleton Laboratory, Chilton Didcot, Oxfordshire, 
0X11 0QX, UK 
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